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Abstract 

We propose a local renormalization group procedure where length scale is changed in spacetime depen- 
dent way. Combining this scheme with an earlier observation that high energy modes in renormalization 
group play the role of dynamical sources for low energy modes at each scale, we provide a prescription to 
derive background independent holographic duals for field theories. From a first principle construction, it 
is shown that the holographic theory dual to a D-dimensional matrix field theory is a (L> + 1) -dimensional 
quantum theory of gravity coupled with matter fields of various spins. The gravitational theory has {D + 1) 
first-class constraints which generate local spacetime transformations in the bulk. The {D + 1) -dimensional 
diffeomorphism invariance is a consequence of the freedom to choose different local RG schemes. 
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I. INTRODUCTION 



Can one prove AdS/CFT correspondence [ll|-|3|]? If so, it will not only give us more insight 
into the precise content of the duality but also open the door to construct holographic duals for 
general quantum field theory. There have been efforts to derive holographic duals directly from 
boundary field theories [|4l-ll6ll. One approach is to build up a bulk spacetime, by introducing 
dynamical sources and their conjugate fields in exchange of decimating high energy modes at 



each step of renormalization group (RG)[18l uM uM- This procedure amounts to an exact change 



of variables from the original D-dimensional fields into the {D + 1) -dimensional variables in 



functional integration, where the length scale becomes the extra dimension in the bulk [| 17 



m. 



The first principle construction provides microscopic justification for the dictionary of the 
AdS/CFT correspondence. If the prescription is applied to a D-dimensional gauge theory, one 
obtains a (D + 1) -dimensional field theory of closed loops which are coupled with a two-form 

n 

gauge field in the bulk lll6ll . The two-form gauge field is an emergent gauge field in the sense that 
its dynamics is solely generated by other loop fields. Because the gauge group is compact, there 
is a topological defect (NS-brane) for the gauge field. Proliferation of the defects describes quan- 
tum tunnelings between different topological sectors. For a sufficiently large A^, the topological 
defects are dynamically suppressed, leading to the deconfinement of the two-form gauge field in 
the bulk. It has been emphasized that those phases that admit 'classical' holographic description 
possess a non-trivial quantum order associated with the spontaneous suppression of the tunneling 
between different topological sectors 016]. This is analogous to the quantum order that is present 
in exotic phases of condensed matter systems with an emergent one-form gauge fieldB, which 
often requires a large number of flavors. 



Despite some progresses, the constructionOSt |l4t |l6|] has an important drawback : it is not back 



ground independent. As a result, it has not been easy to see the emergence of gravitational theory 
in the holographic description. In this paper, we provide a prescription to construct holographic 
duals in a background independent manner. Using the prescription, we show that a /^-dimensional 
matrix field theory can be mapped into a {D + 1) -dimensional quantum gravity coupled with mat- 
ter fields of various spins. This allows one to identify the boundary field theory as a quantum 
theory of gravity in the bulk. 

Here is an outline of the paper. In Sec. II, we start by defining a concrete matrix field theory 
whose holographic dual will be constructed in the remaining of the paper. A theory is defined by 
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specifying sources for all operators allowed by symmetry. A local operator is constructed from 
traces of the fundamental matrix field and its derivatives. Although the field theory is nominally 
defined on the flat D-dimensional spacetime, one can view the theory with spacetime dependent 
sources as a theory defined on a curved background spacetime. In Sec. Ill, we eliminate multi-trace 
operators by introducing a dynamical source and its conjugate field for each single-trace operator, 
where the conjugate fields represent the operators themselves. In particular, the D-dimensional 
metric and its conjugate field that represents the energy-momentum tensor become dynamical. In 
Sec. IV, a local coarse graining is performed where the length scale is increased at a rate that 
depends on spacetime. As the high energy modes are integrated out, non-trivial actions are gen- 
erated for the dynamical sources and the conjugate fields. One contribution is a source dependent 
determinant for the quadratic action of the high energy mode that is integrated out. This includes 
the D-dimensional curvature term of the dynamical metric generated a la Sakharov's induced 
gravity fl20ll . The other contribution represents double-trace operators generated from quantum cor- 
rection, which become a quadratic action for the conjugate field as the double-trace operators are 
removed by another set of auxiliary fields. In Sec. V, we take the advantage of the D-dimensional 
diffeomorphism invariance, which is ensured by the fact that the D-dimensional metric is fully 
dynamical, to shift the D-dimensional coordinates of the low energy modes relative to the coordi- 
nates of the high energy field. In Sec. VI, it is shown that one can construct a {D + 1) -dimensional 
theory as the coarse graining procedure is repeatedly applied to the low energy mode. In Sec. VII, 
we show that the theory in the bulk takes the form of a (D + 1) -dimensional canonical quantum 
gravity if one interprets the extra dimension associated with the scale as a time. In particular, there 
are (£) + 1) local constraints which originate from the fact that the partition function is independent 
of the local RG scheme : the partition function is invariant under the changes in the speed of local 
coarse graining and the D-dimensional shift. From this, it can be shown that those constraints are 
first-class, which generate (D + 1) -dimensional local spacetime transformations. In Sec. VIII, we 
apply the prescription to a simple toy model (0-dimensional matrix model) to illustrate the main 
idea in the simplest setting. In Sec. IX, the difference between the present holographic description 
and the conventional RG is contrasted. In particular, we emphasize the fact that the beta function 
is promoted to a 'Heiserberg' equation for quantum operators in the holographic description. 
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II. MODEL 



A. Matrix field theory 

Consider a matrix quantum field theory defined on a D-dimensional flat spacetime. To be 
concrete, we consider a theory of x real traceless symmetric matrix field $(a;) with the 
global 0{N) symmetry under which the matrix field transforms as an adjoint field. The 'partition 
function' is 



Z[J] = J exp iN^ J rf^x (-^™0^ + V[0™; J 



(1) 



Here O^'s denote single-trace operators constructed from $ and its derivatives. In general, one 
can take {Om} to be a complete set of primary single-trace operators. Here we use the basis where 



Ojn takes the form of 



(2) 



where g + 1 is the order in the matrix field, and {yU^} specifies the spacetime indices. General 
single-trace operators can be written as linear combinations of these operators and their derivatives. 
For simplicity, we assume that there is no boundary in spacetime. Any operator that has overall 
derivatives is removed by integration by part in Eq. ([T}. Throughout the paper, we will use the 
compressed label, say m to denote the full indices, [q, {yU^}] of a single-trace operator. Explicit 
indices will be used only when it is needed. J"^{x) is the spacetime dependent sources for the 
corresponding operator Om- The information on the signature of the background metric is solely 
encoded in the sources. We assume that the spacetime has the Minkowskian metric with the 
signature (—1, 1, 1, .., 1) for with = 0, 1, .., — 1). V represents a multi-trace deformation, 

oo 

V[0^; = 5^ J{">>'H}0^, (9,i..9,.^0„,) {d,..AiOm,) ... [d,..A^,Om,^. 

q=l 

(3) 

where J7^"^*J^'^'^j'^(x)'s are sources for multi-trace operators. All repeated indices are summed over. 

To make sense of the partition function, the theory should be regularized. Here we use the 
Pauli-Villar regularization. Namely, the sources for high derivative terms are turned on in the 
quadratic action for the matrix field to suppress UV divergence in loop integrals. For example, 
one can use a regularized kinetic term, — tr[$ne"«^$], where □ = d/j^d'^. The mass scale M 



in the higher derivative terms plays the role of a UV cut-off. It is noted that the divergence in 
the determinant of the quadratic action is not regularized by the higher derivative terms. In this 
sense, the partition function itself is not well defined. What is well defined is the ratio between two 
partition functions with two different sets of sources where the divergences from the determinants 
cancel. For example, the divergence in the determinant is canceled in correlation functions of local 
operators. 



B. From flat to curved background spacetimes 

Suppose the manifold is endowed with a background metric G^i,. One can define covariant 
operators that transform as tensor density of weight one under coordinate transformations. 



(4) 

where \/\G\ = a/| det Gapl and is the covariant derivative associated with the background 
metric. Any operator Om defined on the flat spacetime can be expressed as a linear combination 
of the covariant operators, 

OUx) = O(G)0^(x), (5) 

where c^{G) is the transformation matrix. It is a function of the metric G^j, and its derivative at 
the position x. For example, 

0[2,H = [0\i,^,u] + ^tu0^2,x]\ , (6) 



where Op.n = ^tr($9^9.<l>), Og^^,, = ^ v^tr($V^V^$), Og,, = j^^H^V^^), and 
r^j, is the Christoffel symbol for the metric G^i,. Therefore, the same Lagrangian in Eq. ([T} can 
be written in terms of these 'covariant operators', 

where 

jG'm^^^ = :r"(x)c„™(G). (8) 

The inverse of the transformation is given by 

OZix) = d^^G)On{x), 

jrr^{x) = J^'"(a;)ci„"^(G), (9) 



with dJ^{G)cJ^{G) = 5^. The multi-trace operators can be also expressed in terms of the covari- 
ant operators and the covariant derivatives of them, 



J 



2 



9=1 



oZ, [v%..v%oZ)i [vZ...v%pZ) ... (v5..v5^o:;^^^) ,(io) 

where j^'^^'^i^^'^"]^ can be similarly expressed as a linear combination of J^^'^'-^^^^^]^ so that Eq. 
(flOl) coincides with Eq. ([3]). The explicit form of the transformation is not important. The factor 
of |G'|"2 is introduced to make the whole expression to have weight one when jG-Amr}A^]} has 
weight zero. 

This means that the original theory defined on the flat spacetime can be viewed as a theory 
defined on a curved spacetime with any background metric. The theory does not depend on the 
background metric because different choices of metric can be compensated by metric dependent 
sources. However, there is a natural choice of metric. We choose the metric G''^°^"^, such that the 
two derivative kinetic term takes the canonical form, that is, = G^'^^^" , where j(o)[2'H 

denotes the source for ^/\G^\tr $ V|^'°' V^*°' $ . This is always possible because one can make 
j(o)[2,MJ^] symmetric in /i and v without loss of generality. In D > 2, there is a unique metric that 
satisfies the canonical condition, = G^'^^^" for a given set of sources (For a proof of this, 

see Appendix A). Such choice of metric is not unique at D = 2 where one needs an extra condition 
to fix the freedom associated with the dilatation. Here we assume that D > 2. In this choice of the 
background metric, the kinetic term takes the canonical form 



C = -N,J\G^\G^^^^'Hr 



+ ..., (11) 



where the same metric is used for the covariant derivative in each tensorial operator and the source 
for the kinetic term with two derivatives. In this sense, an action on a flat spacetime with space- 
time dependent sources defines a natural curved background spacetime. Physically, this amounts 
to measuring the distance on the manifold based on the cost of the action in the limit that the 
amplitude of field is small and the field changes slowly in spacetime. For example, one can set 
the distance between two points to be 1 when the quadratic action that is needed to twist fields 
between the two points is per unit twist and per unit square modulus for each field. In this 
choice of metric, we denote 

0^^\x) ^ rfJ(G(°))0.(x) (12) 



to write 

Here the background geometry is in general curved, but it is non-dynamical. 



(13) 



III. AUXILIARY FIELDS AND GAUGE FIXING 



In the holographic construction[l8l [141, 116|], sources become dynamical in the bulk. Therefore 
one needs to introduce a dynamical field in the bulk for each independent operator. Since multi- 
trace operators can be written as products of single-trace operators, it is convenient to remove the 
multi-trace operators at the expense of making the sources for the single trace operators dynamical. 
We introduce a pair of auxiliary fields for every single-trace operator llj], 

Z = j Dj^^^'^Dp^^^D^ e^^'^^^'^S (14) 

where 

(15) 

Here and pn'^ (x) are D-dimensional auxiliary fields which in general carry D-dimensional 

spacetime indices, g^^ is an arbitrary D-dimensional metric that is used to define a new set of ten- 
sorial operators O^. It is noted that g is in general different from G^^^ The matrix that transforms 
fields defined with different background metrics is given by 

f^''{G^'\g) = d^\G^'^)c-{g), (16) 

which satisfies f{g, g) = /, f {g^^\ g'^'^) f {g^^\ g^'^) = f{g^'\g^'^) and f {g^'\ g^^y) f {g^'\ g^'^) = 
I. j(i)"'s are dynamical sources for single-trace operators, and pn^'s represent the operators 



themselves [Il6ll . In Eq. (fT5l) . j*^^^" is a Lagrangian multiplier that enforces the constraint pn^ = O^. 
Once and pn'^ are integrated out, Eq. (fTST ) becomes Eq. (fUt . 

Since the partition function is independent of the metric g^i,, we can formally integrate over 
different choices of g^y and divide by the volume of the space of metric (gauge volume). 



The resulting theory has a gauge symmetry associated with different choices of the metric g^i,. It is 
emphasized that this gauge redundancy is different from the coordinate redundancy generated by 
the £)-dimensional diffeomorphism. Rather it is associated with different choices of background 
metric in a fixed coordinate system. Under the gauge transformation, the fields transform as 



J 



(l)m 
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9i^u + 



01 ^ 0^^'^ = f-{g + 5g,g)0l 



(18) 



Again we fix the gauge by requiring that the quadratic kinetic term has the canonical form. 
Namely, we will choose the gauge where j(^)[2'H' — g,A"^', Here we have to take into account 
a non-trivial determinant in gauge fixing because the sources are dynamical unlike the case in Sec. 
11. The determinant associated with the gauge fixing can be obtained from the standard Fadeev- 
Popov method. We first define A(j, ^f) such that 



j D5g 5(/7j2'H(^,^ + 5^) - {g + 5gr)A{j,g) = 1. 



This identity is inserted into the partition function, 

DgDj^^^''Dp^n^D<^DSg 



gauge 



By changing the variables, 

jii)^ = 3^'^V^"'ig,g + 5g), 
Pin^ = f^{g + Sg,g)p^:\ 

and by using the gauge invariance of Ci and the facts that 



dj 



(1) 



dPW 
dg 



Aij('^-f:^{g,g + 5g),g + 5g) = A(j«-,^), 



(19) 



Hj^'\g)S {j^'^''fj'''''\g, g + Sg)-(g + 5gr) e^^'°- (20) 



(21) 



(22) 
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we obtain 

Z = 



DG^^WjW'^DPn D^DSg j{i)^ ^(i))^ ^j{i) [2,h _ ^^^d°x £;[J(l),p{l),G{l)M^^,] 

J DJ^'^^DPi'^D^ ^(j(i))g./d«./:;[jW,p(i),j(i)P-],*]^ (23) 



where 



(24) 

Here /^"(0, 1) = G*^^^) and Om^ refers to covariant operators written in the background 

metric G^^'^^"" = j(i)[2.M'^] which is also the source for the kinetic term with two derivatives. In the 
following, we will use G^'^ and jI^-a"^] interchangeably. The determinant becomes 

A{J) = A(J,JP,H) 

D5G 5 (J"/J2'^'^1(G, G + 5G)-iG + SGY") 



DSG 6 



dyr{x) 



5G^^{y) 



SG'^^iy) - bG^"- 



Qliv — j\1,lJ.u\ 



det 



wiicic 5G'f"'(y) — <5G'"'(j/) 



G'=G 



and we used the fact that ^ 



g'=g 



(25) 



_ _ _ _ g'=g 

'^M) '•^^ Kronecker delta function for symmetrized indices, and 6{x — y) is the -dimensional 

Dirac delta function. To obtain the expression in the third line from the second line in Eq. (l25l) . we 

use the fact that there is one and only one solution for the gauge fixing condition (see Appendix 

A). 



IV. COARSE GRAINING 



Now we perform a coarse graining by integrating out high energy modes of the matrix field 
$. Although the sources are also dynamical fields, one can treat them as background fields 
when one integrates out high energy modes of $. We focus on the functional integration of the 
original dynamical field $ which is coupled to the sources J^^^'^, 

Z^ljW] = j /)$e*/'i°-Kr{4.A^.(i)*)+^.(i)I*]], (26) 

9 



where J^jw is the kernel for the quadratic action that includes the two and higher derivative 
terms, 



n=3 

and Uj(i) [$] includes all other single-trace operators, which are at least cubic in There is no 
operator linear in $ because $ is traceless. The sources for the higher derivative terms in Eq. (|27]) 
have the engineering scaling dimension [jI^'A^i - •A'nlj _ _ 2)^ and the mass scales associated 
with the sources play the role of UV cut-offs. It is interesting to note that there are in general 
many scales. Moreover, the cut-off scales are fluctuating because the sources are dynamical. Each 
configuration of jP'A'i ^"1 describes a theory of $ with a different set of UV cut-off scales. We 



perform a real space RG transformation [|2 ll 12211 by lowering some of these energy scales. For this, 
an auxiliary traceless real symmetric matrix field $ is added to the original theory, 

where Ad is an arbitrary kernel for the auxiliary field. Here (^+^K^~i) the number of indepen- 
dent components of a real traceless symmetric matrix. We go into a new basis (p and (p, 

<l>(x) = 0(x) + 0(x), 

<l(x) = Jdy (A(x,y)0(y) + i?(x,y)0(y)), (29) 

where the functions A and B are uniquely chosen from the conditions that the low energy field (j) 
and the high energy field (p do not mix at the quadratic level, and that the low energy field has a 
set of UV-cut off scales which are smaller than those for the original field $. Then the partition 
function takes the form of 

where the quadratic action for the low energy mode is given by 



n=3 

Here the rescaled sources are 
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where c„ is a set of constants which determine how we rescale the set of UV cut-off scales. Chang- 
ing the sources for high derivative terms of (p in this way is equivalent to lowering the UV cut-off 
scale associated with the source j(i)[2.Mi, -,Mn] ^ factor of e"'^""'^'''^, where a^^^ is the rate at 
which the UV cut-off is lowered, and dz is an infinitesimal constant. It is noted that one can 
choose different speeds of coarse graining at different points in spacetime, and a^^\x) is in gen- 
eral position dependent. This is a local RG procedure where the speed of coarse graining is space- 
time dependent. Specifying {cn} corresponds to choosing a particular RG scheme. One natural 
choice would be Cn = (n — 2) which reflects the fact that •a'i] has dimension — (n — 2). 

This amounts to lowering all UV cut-off scales in the same way. However, this choice is not 
ideal for our purpose : the ratio between the determinants of the original field and the low energy 
field, [det J^ jay det is divergent. To illustrate this, let us consider the simple case where 

j(i)[2,/ii...Atn] a^'^\x) are independent of x. In this case, the ratio between the determinants is 
given by 

/Y^ r 7" 7'{l)[2,Ati---Mn]i. k 
V I~ja)[2........i/'""A^^'"' - (33) 

which is divergent if c„ = (n — 2). In the conventional RG procedure, the determinants do not 
play an important role, and one can ignore the divergent determinant in computing beta functions. 
In our case, the determinants are important because they provides a non-trivial action for the 
dynamical sources. This difference comes from the fact that sources are dynamical fields in our 
approach, instead of constants. In order to avoid the UV divergence, we choose, among many 
other choices, the following prescription, 

Cn = [n — 2) for n < ric, 

= for n > ric, (34) 

where ric is a large but fixed number. For sufficiently large momenta, Eq. (l33l) becomes 

d'^k . ,,,,,, 7- , (35) 

which is finite. This is a well defined coarse graining procedure, where we are rescaling the sources 
for higher derivative terms upto the n^-th order to eliminate high energy mode while avoiding the 
divergence in the ratio of the determinants. In a sense, we are performing a coarse graining with 
two sets of scales. The first set of scales associated with the high derivative terms up to the nc-th 
order plays the role of the usual UV cut-off that is rescaled to thin out high energy modes. The 
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second set of scales associated with the high derivative terms with more than Uc derivatives cuts 
off the UV divergence in the ratio of the determinants. The conventional scheme is reproduced 
when Uc is taken to be infinite. It is emphasized that the specific form of rescaling in Eq. (|34l ) is 
not important. There exist many other schemes that regularize the divergences in the determinants. 
What follows below is independent of the specific choice. The propagator of the high energy mode 
is given by the difference between the propagators of the original field and the low energy field 

= M-l,-M-l,y. (36) 

Therefore the propagator of the high energy mode is 0{dz). 




FIG. 1: Two ways of generating quantum corrections to the Linear order in dz. Each circle denotes trace 
of a chain of matrix fields. Solid lines represent chains of low energy fields and each dashed line represents 
a high energy field, (a) Contraction of a pair of high energy fields within a single-trace operator generates 
two singe-trace operators (the first and the third) and one double-trace operator (the second). In the large 
A'^ limit, only the second term is 0{N'^). (b) At the quadratic order, one can fuse two single-trace operators 
each of which contains one high energy mode. This leads to one double-trace operator and two single-trace 
operators, all of which are 0{N'^). 

Integrating out the high energy mode, we obtain an effective theory for the low energy mode, 

(37) 

where 5J^^^ is the quantum correction to the sources for the single-trace operators, and 
denotes double-trace operators generated from quantum corrections. Because the propagator of the 
high energy mode cj) is of the order of dz, only two diagrams contribute to the quantum corrections 
to the linear order in dz. The first contribution comes from contracting two high energy fields 
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within one single-trace operator as is shown in Fig. [11(a), 



5C = -ArJ|G«|j(^) 



tr 



(V{.^}0) .. (V 



V 



..(V 



V{M^}0) (V 



..(V 



(38) 



whei^ (V|,.>0) ^ (Vf Vg^.Vg;V) and 

j D(j)0 e^^/'^''^tr(^'5<^) 
^ *^ J L)0 ^■^^('^'50) 

Using the expression for the propagator of the high energy mode, 



(39) 



< (j)ab{x)(j)cd{y) >0 = i 

one obtains 

5C = 



8N 



i^ac^bd + ^ad^bc) — J^^ab^cd 



(40) 



tr 



0(V|^i}0)..(V{^.-i|0) (V 



(V 



tr 



..(V 



.(41) 



Here we drop the contributions that are sub-leading in Integrating by part, if necessary, one 
obtains operators of the form, 

jV2^W^W ^^(1) ^(i)^^^(l) 



-(vSo«)(v;:.\oW) 



p,9=0 



5 ,i)jW"^"-'t^i'--''^p''^i'--''^i 



(v«..v«o(i))(v«..v«o«). 



(42) 



This expression represents double-trace operators when both Om and On ^ are non-trivial. In the 
special case where Om = 1 with p = (or On^ = 1 with q = 0), it reduces to single-trace 
operators. The special case occurs when two high energy fields are adjacent to each other in a 
chain of matrix fields within a single-trace operator, as in jI'^'H ^/[(^(i) | j^tr[(f)(f)(j)V ^(j)] . The 
second contribution comes from fusing two single-trace operators as is shown in Fig. \T\(h), 

tr 



(43) 



tr 



(f) (V{^i}0) .. (V{^.-i}0) (^V{^.}0j (V{^>+i}0) .. (V{^<,}^ 
(V|,i}0) .. (V|,,-i}0) (V|,,|0) (V|,,+i}0) .. (V{,p}0] 



y/ (f) 
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Contracting the high energy modes, one obtains both single-trace and double-trace operators, 

5£' = / d^'y y/|G'(i)(x)|^|G(i)(y)|jW[^+i'^^5}l(a;)J«[^+i-H}](y) x 



Vl^.^V\^4S-\x,y) + S-\y,x)]{-2ir\<p (V{^i}0) .. (V{^.-i|0) (V{^«+i}0) .. (V 



X 



tr 



+A^tr 
-fA^tr 



V{^1}0 j .. (^V{^.-1}0 j (^V{^j + l}0 j .. (^V{i.P}0 

(V{^i}0) .. (V{^.-i}0) (V{i..-i}0') (^V{^j-2|0'j .. (^V{^.+i}0'j (V{^.+i}0) .. (V{/,.}( 
(V{^i}0) .. (V|^«-i}0) ( V{^,+i}0 j ( V{^,+2|0') .. ( V|^,-i}0 j (V|^«+i}0) .. (V{^<,}< 



(4 



In this expression, it is understood that (j) = (j){x) and (j) = (f){y). Although the quantum correc- 
tions appear to be non-local, the propagator for the high energy mode decays exponentially in real 
space, allowing one to do a gradient expansion. The scale that controls the gradient expansion is 
the UV cut-off which is set by the dynamical sources. This results in local double-trace opera- 
tors of the form in Eq. (|42] ) and quantum corrections to single-trace operators. The determinant 
[det A^j(i)/ det ' 4 - gives rise to a Casimir energy that depends on the source J^^\ 

The Casimir energy provides a 'potential' energy for the dynamical sources while the double- 
trace operators becomes a quadratic 'kinetic' term for the conjugate fields as will be shown in Sec. 
VI and VII. 

Note that even though the action for $ in Eq. (|24l) has only single-trace operators, double-trace 
operators are generated in the renormalized action [liol 111! . Triple or higher trace operators are at 
least order of dz"^ and can be ignored in the small dz limit. After exponentiating the determinant 
into a quantum action for the dynamical sources, the total Lagrangian for the low energy field and 
dynamical sources can be written in the following form, 

£2 = iV2|v[/^"(0,l)P«;:r(°)'^'"'^'H}] 



(Vi;i0(i))(V«0i^))|. (45) 



+ (j(i)- _ j(o)-/„"(0, 1)) P(i) + 5,(.)£[J(^)^ 



5^(1) j(i)™«M{^} 



Here /rf^x = -i ^^^^j^^"^^ \n[det Mjw detMj',^]. We use the same notation 

Om to represent single-trace operators constructed from the low energy field 0. d^w jW^^t^^ 
includes both the change in the quadratic sources caused by lowering UV cut-off in Eq. (|32] ) and 
the quantum corrections to the single-trace sources in Eqs. (|4TI) and (l44l) . denotes 
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the sources for double-trace operators generated from quantum corrections in Eqs. (HTT ) and (|44] ). 
In the Casimir energy and the quantum corrections to the sources, it is enough to keep only those 
contributions to the order of dz. Therefore, 6^wC[J^^^% J^^)™^^^ and J^^)™"^'^^^} are 
linear in a^^^dz. In general, all terms that respect the D-dimensional diffeomorphism invariance 
are allowed in the Casimir energy and the quantum corrections to the sources, 

= d^aW(a;)^^G(l^{Co[J«] + Cl[J«]7^ + ...}, (46) 
5^(1) J(')'"^'^^ = dza^'\x)A'^^''^[J^% (47) 
= rf2a(^Ha;)5""^'^>^'^>[j(^)]. (48) 

Here is the Casimir energy which includes the cosmological constant Co and the 

LJ-dimensional Ricci scalar TZ. It also includes higher order terms in the curvature tensor and 
derivative action for other sources. Since the Casimir energy comes from the deter- 

minant of the quadratic operator of $, it depends only on the sources for the quadratic operators. 
On the other hand, in Eqs. (|47l) and (|48T ). cubic or higher order operators are also renormalized 
by quantum corrections. Cq, Ci, A"^^^^ and B'^'^^^^^'^^ are finite functions of the sources 
and their covariant derivatives. Note that one needs to include descendant operators for quantum 
corrections in Eq. (|45T ) to express 6^(i)J^^^'^^^^ and S^(i)J^^^'^"^^^^''^ as linear functions of a^^'^ 
without derivative as in Eqs. (|47] ) and (|48] ). 

n 

It is noted that the partition function is independent of a(^)(x)[8]. This is because is an 

arbitrary function introduced to change the length scale in the coarse graining procedure. One can 
choose any speed of RG without affecting the partition function. If one modifies the parameter a^^^ 
to a^^^ + Sa, the quantum corrections will be modified accordingly, exactly undoing the changes 
caused by 5a. In other words, one has to add quantum corrections (counter terms) so that the 
partition function computed from the low energy effective theory with a lower UV cut-off is equal 
to the one computed from the bare theory with the original cut-off. This is nothing but the cut-off 
independence of the partition function in the Wilsonian RG. The fact that the partition function is 
independent of the choice of a^^^ (x) will become important later to obtain the {D + 1) -dimensional 
diffeomorphism invariance in the holographic description as will be discussed in Sec. VII. 
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N^dz 
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gh energy 



FIG. 2: The D-dimensional coordinate of the low energy mode is shifted infinitesimally relative to the 
coordinate of the field defined at high energy. 

V. SHIFT 



One key difference of the present construction from the conventional RG procedure [|2li |22|] is 
that the sources (and the conjugate fields) that are coupled with the low energy field are dynamical. 
In particular, the low energy field is covariantly coupled with the fully dynamical D-dimensional 
metric. Therefore, there is a freedom to choose the coordinate system for the low energy mode 
without modifying the form of the Lagrangian. To take advantage of this extra freedom, we will 
choose a new coordinate system for the low energy mode which is infinitesimally shifted along a 
D-dimensional direction compared to the original spacetime [|l2n . The point is that one does not 
have to use the same D-dimensional coordinate for the low energy field as the one for the high 
energy field as is illustrated in Fig. |2l For this, we single out the action for the low energy mode. 



S, = d^x 



X T(l)mn{u}{i/} 



(49) 



and change the variable, 

4>{y) = 0(x) (50) 

with = + N^^'^^{x)dz. In the new coordinate system, the metric and the sources are trans- 
formed as tensors. The operators transform as tensor densities of weight one. 
To the linear order of dz, we have 



Si = I d^y 



_ 7(l)m/)(l) _ r T(l)mMv(^) + 



(51) 
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where Om's are the covariant operators constructed with the covariant derivative V^^^ with the 
new metric 

= (52) 
These operators are coupled with the new sources given by 



J(l)[<?'{''5}](y) = 



n 



J(1)[<?,H}](^)_ (53) 



We can ignore the shift in 5^(1) and ^^.(i) because they are already order of dz. 

The Lagrangian in the shifted coordinate becomes 

+ (vaQi^))(vaoi^))}, (54) 

where Sj^wG^^^"" = G'^^)'*'^ - G^^^'"' and 

As was the case for a'^^\ the partition function is clearly independent of N^^^^^, because different 
choices of shift merely corresponds to different choices of coordinate system for the low energy 
mode. This completes one cycle of the RG procedure. We have a theory of the low energy field 
coupled with dynamical sources whose fluctuations are controlled by the action generated from 
the high energy mode. 



£2 = 7V2|v-[/^-(0,l)Pii); J(°)'^™'>'H>] 
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VI. CONSTRUCTION OF BULK THEORY 

Now we repeat the procedures in Sees. III-V. Another set of auxiliary fields are introduced to 
remove the double-trace operators for the low energy fields followed by the gauge fixing to obtain 



where we use the notation $ for the low energy mode to avoid introducing a new notation for the 
low energy mode at each step of RG, and the Lagrangian is given by 



+ (J(l)n _ j(0)-^^n^Q^ p(l) ^ 5^^,,£[J(1)-] 



Here 0„ s represent covariant operators constructed with the metric G^'^^^'^ = j(2)[2,/xi/] 
in Eq. (l57l) is the Jacobian generated from the gauge fixing. As was done in Sec. IV, high en- 
ergy modes are integrated out with a spacetime dependent coarse graining rate a^'^\x) to generate 
the Casimir energy for J^^)" and another set of quantum corrections to single-trace operators and 
double-trace operators, which are proportional to a^'^\x)dz. This is followed by another infinites- 
imal shift along the D-dimensional direction N^'^^^{x) as in Sec. V. Note that the a^'^\x) and 
N^'^^^(x) are independent of a^^\x) and N^^^^{x). Namely, we can choose different rate of coarse 
graining and different shift at each scale. The double trace operators are again removed by intro- 
ducing a third set of auxiliary fields j(3)m, p(?) 

If we repeat these steps L times, the partition function can be written as 



(57) 




(58) 




(59) 
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where 



£4 = 7W[/^"(0,1)PW; J(°)'^'"'>'H>] 

L 

+N' ( jC+i)" - / + 1)) Pi'+^) } 

L 
1=1 

_^2j(L+l)nQ(L+l) (gQ) 

Here it is understood that J(°)" = 

Now we first take the Umit with dz ^ and L ^ 00 with fixed zl — Ldz, where z = Idz 
becomes a continuous coordinate that labels the length scale in the range < z < zl- J^^^^{x), 
Pn\x), a^''\x) and N^^^'^{x) become D + 1-dimensional fields J^{x,z), P^{x,z), a{x,z) and 
N^{x, z), respectively. Then, we take the — > 00 limit, which amounts to taking the low energy 
limit where one push the RG procedure to the IR limit. In this limit, the partition function becomes 



Z[J] = J VJ{x,z)VP{x,z)M{J) e 



SuvIP{x,0)]+SIJ{x,z),P{x,z)]+Sir[J{x 



J{x,0)=J{x) 

(61) 



where 

00 

VJ{x,z)VP{x,z) = Y[[DJ^^\x)DP^^\x)] , 

1=1 

00 

M{j) ^ n^('^^'^ 
1=1 

Suv = J d^'x V[P^(a;,0); J(°)'{"^^>'H>], 

S,r, = -ilnj D* e-"''/''°'-''"<'.~)o2-~'. (62) 

Here Suv and Sjr are the actions defined at the UV (z — 0) and the IR (z — oo) boundaries 
respectively. 5" is the bulk action. J" = (9^ J" + J'^A^) is a 'covariant derivative' for a vector 
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J" defined in the space of operators with the connection, 



A::,{x,z) = dJ^{G{x,z),G{x,z')) 



'5Gt^-{yy ^^^^ 

Note that is not related to the covariant derivative in the D-dimensional spacetime V^. H and 
T-i^ are given by 

\/\G\ 

(64) 



G\\Co[J{x)] + Gi[J{x)]n+ ... 
1 



2^/\G\\/' 

E 

[g,{/.n]7^[2,H 



p 



X 



5f 



6G^"'{x) 



1 



a, 6 



j[g,{M}/^2-M6-i'^M6+i---}] p 



[q,{nll^il...fi'^_j^fj,Hl }] 



+ (V^J[^'^^5>1)P[,, 



(65) 



with /i = 0, 1, ..,(£> -1). 

Starting from the D-dimensional matrix field theory, we obtained a {D + 1) -dimensional theory 
for dynamical sources and operators. The extra dimension parameterized by z represents the length 
scale in the RG. One can choose different speed a{x, z) of RG at different points in spacetime and 
scale. Therefore z is not a gauge invariant quantity. What is gauge invariant is the length scale 
whose infinitesimal increment is given by dr = adz. In this sense, the physical length scale can 
be viewed as a proper length along the direction of the extra dimension. This is illustrated in Fig. 
[3l The theory in the bulk has the dynamical metric and its conjugate momentum as dynamical 
degrees of freedom. Therefore it is natural to expect that the bulk theory is a gravitational theory. 
In the next section, we will see that the theory in the bulk indeed respects the {D + 1) -dimensional 
diffeomorphism invariance. 



VII. HAMILTONIAN GRAVITY 

If the RG scale z is interpreted as 'time', the theory can be viewed as a Hamiltonian system 
where the sources J™'s play the role of 'coordinates' and the operators P^'s are the 'momenta'. 
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(a) 



X 



(b) 



X 



FIG. 3: (a) Bulk spacetime made of the D-dimensional boundary spacetime and the semi-infinite line that 
represents the length scale in the RG procedure. Each step of coarse graining, say the Z-th step, generates 



a set of D-dimensional fields ( j(')"(x),Pi \x)] that represent dynamical sources and operators at that 



scale. These fields are combined into {D + 1) -dimensional fields ( J"(x, z), z)) in the bulk, where 
the extra coordinate is given by z = Idz. Each 'vertical' line traces the positions of the bulk fields which 
are generated from the original field variable ^{x) at each x in the boundary spacetime. The spacetime 
dependent shift N^{x, z) causes the bulk fields to have different D-dimensional coordinates from that of 
$(2;). Each 'horizontal' line represents the manifold in the bulk spacetime with an equal z coordinate. 
Because the speeds of coarse graining are in general different at different points in spacetime, two points 
within the manifold with an equal z do not in general have the same proper length along the extra dimension, 
where the proper length is the scale in the RG. (b) The same bulk spacetime where the coordinate z is used 
instead of the proper length along the extra dimension. The vertical lines have the same meaning as in (a). 
Each horizontal line represents the manifold with an equal proper length, that is, the set of points with the 
same length scale in RG. Note that an horizontal line that is concave upward in (a) is concave downward in 



The sources and operators are conjugate to each other as expected. The 'Hamiltonian' is given by 



where N^{x, z) = a{x, z) and = Ti. Note that the 'time' x^ = z is different from the real 
time x° in the boundary field theory. The Hamiltonian in Eq. (|66l ) generates the evolution along the 
time x^ associated with increasing length scale of the system, not along the real time In this 
sense, one can regard the Hamiltonian as a generator for a quantum beta function. One difference 



(b). 




M=0 



(66) 
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from the usual Hamiltonian system is that the 'covariant' derivative Vz is used in the action. The 
non-trivial connection originates from the fact the fields J"{x,z), Pn{x,z) defined at different 
length scales in general have different metric. Because the definition of the covariant operators 
and their sources is tied with the metric, a change in metric effectively induces changes in all 
sources. Physically, the momentum canonically conjugate to the metric is the energy-momentum 
tensor given by 7r[2,^i.] 
tensor besides P[2,^iu] = 



There are many other contributions to the energy momentum 

Id /\G\ 
^tr[$V;,V,.$] because metric enters not only in ^^*^tr[$V/,V,.$] but 



also in the definition of all other covariant operators. This suggests that P[2,^v] is not the canonical 
momentum of the metric, which is also reflected in the non-trivial measure A (J) in the functional 
integration, and the unconventional form of the shift for the metric in Eq. (|65T ). In order to go 
to the canonical basis, we define a new momentum for the metric and keep the same conjugate 
momenta for all other variables. 



(67) 
(68) 



Pm, form 7^ [2,iJ,u]. 

The last term in Eq. (|67T ) takes into account the metric dependence in general operators. The 
Jacobian from the change of variable. 



5P[2,f,u]iy) 



[2,a/3] 



det 
A(J) 



t''^)Kx - y) + ny) 



SfJ''"'\y) 
5G'''^(x) 



-1 



-1 



(69) 



exactly cancels A (J) in the measure. The partition function and the action takes the canonical 
form[23] in the new variables. 



Z[J] 



VJ{x,z)V7r{x,z) e 



■fi)]+S[J{x,z),ir{x,z)]+SiR[J(x,oo)] 



J{xfi)=J(x) 



where 



S = d^xdz 



(9^ J")7r„ - a{x, z)'H - iV^(a;, z)n^ 



(70) 



(71) 



Note that is replaced by the usual derivative in the canonical variables. Moreover, the 'mo- 
mentum constraint' T-L^ that generates the D-dimensional shift takes the standard form. 



(72) 
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It is noted that J^^'s and vTm's are D-dimensional contra-variant tensors with weight zero and 
covariant tensor density with weight one, respectively. To obtain the 'Hamiltonian constraint' T-L, 
one has to convert Eq. (|67] ) to express P[2,f^u] as a linear combination of n^s and plug in the 
expression to Eq. (|64]) . Since the full expression is complicated, we focus on the metric and its 
conjugate momentum. Among many other terms, "H includes the linear and quadratic terms for 
the conjugate momentum, the cosmological constant and the D-dimensional curvature, 

H = A^'^[J(a;)]7r[2,H ^^p^7r[2,H^[2,A.] 

V 1^ I 

-^\[Co[J{x)] + Cl[J(x)]7^} + (73) 

where ... represents the higher dimensional terms that involve covariant derivatives of vr and the 
curvature. Cubic or higher order terms in 7^[2,fj.u] are not allowed because at most double-trace 
operators are generated out of single-trace operators at each step of RG. The linear term in the 
conjugate momentum arises because the operators that are quartic in $, such as -^tr[$^V^V;y$], 
renormalizes the metric through the quantum correction in Eq. (|4T1) . It is interesting to note that 
the kinetic term for the conjugate momentum originates from the beta function under the RG, 
while the potential term for the metric originates from the Casimir energy. Besides the dynamical 
gravitational mode, the theory also includes other degrees of freedom, including the higher spin 
fields for j^tr[^V /j^^.-.V ^^^] and the fields associated with the single-trace operators that are 
cubic or higher order in $ . As was noted in Sec. IV, the latter fields do not have the bare 'potential 
energy' because the Casimir energy is independent of those fields. However, they do have the 
quadratic kinetic term in general because double-trace terms are generated for those operators 
under the the RG. Although the bare action for those higher order sources are ultra-local along the 
D-dimensional space, potential terms that involve derivatives along the D-dimensional space will 



be generated dynamically, as other heavier fields are integrated out in the bulk [|l6ll 



In the large limit, the bulk fields become classical. In particular, non-perturbative fluctuations 



of the bulk fields are dynamically suppressed fll6ll . The on-shell action in the bulk computes the 



partition function of the original matrix field theory in the large limit. The classical equation of 
motion is given by 

a,J" = {J",H}, a,7r„ = {7r„,H}, (74) 
where the Poisson bracket is defined by 



{A,B} = I d^x 



5A 6B 6A 6B 



(75) 
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To solve the equation of motion, we need another set of boundary conditions besides J"(x, 0) = 
J7"(a;). The second set of boundary condition is dynamically imposed by the regularity condition 
in the IR limitLSJ. At the saddle point, we have 

Z[J{x)] = e^(^^^+^)Z[J(x, oo)], (76) 

where the bulk action is evaluated at the saddle point configuration, and J{x, oo) is determined 
from the condition that the bulk action is finite in the IR limit. At the first glance, this expression 
does not seem meaningful because both Z[J'{x)] and Z[J(x, oo)] are not well defined due to the 
divergent determinants. However, correlation functions are well defined because the divergences 
in the determinants cancel. Correlation functions of local operators are given by 

z" d''Z[j'{x)] 



(Xi ) 0„2 (0:2) . . . On^ {Xn^ ) ) 



where 



(77) 

w=0 



j'-(x) = J"(x) + J2^'''SlS{x - Xi). (78) 



If non-local sources are turned on, both the saddle point solution in the bulk and J(x, 00) are 
modified. For perturbations that are localized both in space and time in /^-dimensions as the one 
in Eq. (fTSl) . it is expected that J{x, 00) does not depend on the perturbation, w"'\ This is because 
local perturbations are always irrelevant and die out in the IR limit in unitary theories. Therefore, 
the contribution from Z[J{x, 00)] will drop out in Eq. <^^}. 

It is noted that the D-dimensional general covariance does not allow mass term for the metric. 
This does not preclude the possibility that metric fluctuations become massive as other fields are 
condensed, breaking the D-dimensional Lorentz symmetry. Here we consider the simple case 
where the sources J"^{x) at the boundary respect the D-dimensional Lorentz symmetry, and the 
vacuum does not break the symmetry spontaneously. In this case, the saddle point configuration 
of the bulk fields will also respect the D-dimensional Lorentz symmetry. Mathematically, this 
means that the only tensor that has a non-zero spin and has a non-zero expectation in the bulk is 
the metric. Then, the tensors for the conjugate momentum in Eq. (1731) take the form. 



A^'^lJix)] = A[J{x)]G^"', 
5^'^^"[J(x)] = Bi[J{x)]G'"'G^'' + B2[J{x)]G^'^G'"' (79) 

at the saddle point, where A[J(x)], Bi[J{x)] and B2[J{x)] are scalar functions of the sources and 
their covariant derivatives. 
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Although the theory in the bulk is a quantum theory of dynamical metric in ( D + 1 ) -dimensional 
space, it is not clear whether this theory has the diffeomorphism invariance in the bulk, which is 
the key property of gravitational theories. In the canonical formalism, the {D + 1) -dimensional 
diffeomorphism invariance would show up as (D + 1) first-class constraints. If A[J(a;)], Bi[J{x)] 
and B2[J{x)] were just constants, they have to satisfy specific conditions in order for the Hamilto- 
nian constraint to be first-class. For generic values of A, Bi and B2, the Hamiltonian constraint H 
is not first-class, in which case the theory does not have the full {D + 1) -dimensional diffeomor- 
phism invariance. Given that the coefficients are dynamically determined, it seems highly unlikely 
that they have the saddle point values of the fixed ratio at all points in the bulk independent of 
JT". However, we have to be more careful here because the present theory is not a pure gravita- 
tional theory. As a result, A, Bi and B2 depend on other dynamical fields which themselves have 
non-trivial Poisson bracket with their own conjugate momenta. Namely, we can not just replace 
A, Bi and B2 with the saddle point values when we determine the nature of the constraint. In 
other words, one should compute the Poisson bracket among the constraints, treating all dynami- 
cal fields on the equal footing. Instead of computing the Poisson bracket explicitly, here we use a 
simple argument to show that all {D + 1) -constraints are first-class. 

As was emphasized in Sees. IV and V, the partition function does not depend on the choice of 
the lapse N^{x, z) = a{x, z) and the shift N^{x, z). From the fact that the partition function is 
independent of N^{x, z), we obtain 

Therefore the lapse and the shift play the role of Lagrangian multipliers which impose the local 
constraints, 

7{ = 0, n^.^Q (81) 



inside the bulk spacetime. Since the above equality holds at any time z, we have 

d_ 

dz 

In order for this to be true for any choices of N^{x, z), we have 



<nM{x,z) >= J d^'y N^' iy,z) {{nMix,z),nj^>{y,z)}) ^0. (82) 



{HMix,z),nM'iy,z)} = (83) 

at the saddle point. This implies that the (D -I- 1) constraints are first-class classically. These 
constraints generate local spacetime transformations in the bulk. The Hamiltonian constraint H 
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generates the transformation. 



x'^^x'', z ^ z + dlN^{x,z), (84) 
whereas the momentum constraint H^i generates 

x^ ^x^ + dW{x, z), z^z. (85) 

A general (D+l)-dimensional diffeomorphism generated by a combination of the two corresponds 
to choosing a different prescription for the local RG procedure. 

VIII. A SIMPLE EXAMPLE 

Because the construction is rather complicated, it will be useful to apply the prescription to a 
simple toy model to illustrate the backbone idea. Here we provide an explicit construction for the 
simplest possible matrix model : 0-dimensional matrix theory. The partition function is given by 



exp 



(86) 



i[-Nj2 ^"tr($") + ArV[tr($") /N] 

. \ n=2 

where $ is a real traceless symmetric matrix and l/[tr($") /N] is a general non-linear function of 
single-trace operators which may be expanded as 

oo 

V^Yl Yl A^~^J"''''"'-''"''tr($"^i)tr($'"2)..tr($"^). (87) 

q=2 mi,m2,...,mq 

Because there is no spacetime, the partition function is given by a single matrix integration. J7"'s 
(^j-mi,m2,..,mg'g-j rcprcscnt the sources for the single-trace (multi-trace) operators. We assume that 
the sources have small imaginary components so that the integration is well defined, e.g., ImJ"^ — 
— e for even n; ImJ"^ — for odd n. 

To remove the multi-trace operators in V, we introduce a pair of auxiliary fields, J^^^'", Pn^^ 
for each single-trace operator, 

Z ^ J dJ^^^'^'dPj^^^d^ e'^\ (88) 

where 

^1 = iV^ { [Pi'^ - ^) - + V[P^^^] } (89) 
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jmi,m2,..,m,pWpW j^ic contours of P^^^'s are along 



the real axis, but the contours of J^^) " are chosen slightly off the real-axis as 
which guarantees that integration for $ is well defined. Now we have only single-trace operators 
for $ which are coupled to the dynamical sources J*^^)'". P^^^ is the conjugate variable which 
corresponds to the single-trace operator, -^tr($"). This can be seen from Eq. (l89l) where 
plays the role of a Lagrangian multiplier which enforces the constraint, P^^^ = -i:tr($"). 

There is only one operator which is quadratic in $. We use its source J*^^)'^ as a scale to generate 
a renormalization group transformation. Because < $^ >~ l/J*-^^'^, we can regard l/J^^^'^ as a 
UV cut-off, and generate RG flow by lowering 1/ J*^^^'^[18|]. The fact that 1/ J*^^)'^ plays the role of 
a UV cut-off can be understood from the observation that with a smaller 1/ J*^^)'^ the fluctuations 
of $^ decreases. Using the method described in Eq. (1291 ). the original matrix field $ can be written 
as a sum of the low energy field (p and the high energy field 0, 



(JV+2)(iV-l) 



JS2 



(90) 



where 



00 

-Arj(i)'2'tr(02) _ Nmhr{4>^) -NJ^J^^^'M 



+ 



(91) 



n=3 



with J(l)'2' = J(l),2e2a(l)d. ^nd m2 



J(l).: 



2a(i)d2" f^^^^ '^^ infinitesimal parameter and a'^^^ is 
the rate at which the UV cut-off is lowered in the first step of RG. Integrating out the high energy 
mode, one obtains the effective action which includes a Casimir energy, quantum corrections to 
the single-trace operators, and double-trace operators. 



^3 = iv2|v[pW] + (j(i).-_j7")pa)| 



(92) 



where 

p"«[j(i)] 



1 



(/ + l)(n + 1) J(l).'+1 _ !(!i±i±^ j(l),n+i+2 



(93) 



2J(i)'2 

In this 0-dimensional matrix model, the Casimir energy is a constant independent of the sources. 
However, in higher dimensions, the Casimir energy is in general a function of dynamical sources. 
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including metric as we saw in Sec. IV. Now we introduce another set of auxiliary fields J^^^'", 
Pk^' to remove the double-trace operators that are generated from quantum corrections. Then (j) is 
again divided into the low energy mode and the high energy mode, by rescaling J*^^^'^ by e^"^^^'*^^. 
Integrating out the high energy mode generates double-trace operators, which are removed by 
another set of auxiliary fields. Repeating these steps, one can write down the original partition 
function in terms of the integration of the dynamical sources J^'')'" and the conjugate variables 
Pn"* introduced at each step of RG. As dz — )■ 0, the discrete RG step becomes a continuous 
dimension, and jC^)'", Pn\ a'^^^ become functions of z : J'^{z), Pn{z), a{z). The original 0- 
dimensional theory in Eq. (|86] ) is now mapped into an one-dimensional theory. 



Z[J] 



Vr{z)VPn{z) e 



i(^Suv[Pm+S[J(z),P{z)]+SiR[J{T)] 



(94) 



where 



Suv = iW[P^(0)], 



S = [ dz [{d,r)Pn - aiz)^ 
Jo 



S 



IR 



—i In 



(95) 



Here T is the RG 'time' at which we stop the coarse graining procedure. This creates a boundary 
dX z = T and a boundary action Sir. The partition function is independent of T, and one can take 
T — > oo to push the TR boundary' to infinity. The 'Hamiltonian' in the bulk is given by 

.(iV + 2)(iV- 1) 



2J2 



2iV2 



2rP2 



n>2 



^ 2N \ N 



,fc,«>3;fe+/=2+n 



2J2 



(96) 



J,n>2 l+n>2 

For P„ with n < 2 we use the convention, Pq = 1 and Pi = 0, which reflect the fact that 
-^tr(/) = 1 and ;^tr(0) = 0. This is the one-dimensional holographic description for the 0- 
dimensional matrix theory. 

From the expression in Eq. (|94] ). one immediately realizes that the partition function can be 
viewed as a transition amplitude of a quantum mechanical system if z is identified as time and H as 
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Hamiltonian. In the Hamiltonian interpretation, the sources and the conjugate fields are promoted 
to quantum operators, and satisfy the commutation relation, 

[J\Pm]^i^,5l, (97) 

where [A,B] = AB — BA is the usual commutator and plays the role of the Planck constant. 
The partition function can be written as 

Z = < mf\e-'^^i^dza(z)H^^,, ^ (98) 
where the initial and the final wavefunctions are given by 

< J"|^^ >* = e^^^«[-^"l. (99) 

Since Hamiltonian is not Hermitian, the evolution is not unitary, which is consistent with the 
fact RG flow is irreversible. Note that a{z) becomes the lapse function along the time direction. 
Moreover, the partition function is independent of the choice of a{z). Choosing a different a{z) 
amounts to using a different parameterization along the RG flow. The independence of the partition 
function under the reparameterizaion of RG flow is nothing but the diffeomorphism invariance of 
the bulk theory. This one-dimensional diffeomorphism invariance is expressed in terms of the 
constraint, 

Sa{z) 

However, Hl^i > does not identically vanish because the boundary at UV explicitly breaks the 



diffeomorphism invariance. 



IX. HOLOGRAPHIC (QUANTUM) RG VS. CONVENTIONAL (CLASSICAL) RG 

Finally, we compare the holographic description with the conventional RG. In the present con- 
struction, the effective action contains only single-trace operators at all energy scales. This greatly 
simplifies the RG procedure. The price one has to pay is that one has to promote the sources 
of the single-trace operators to dynamical variables. In other words, the couplings are not mere 
constants any more, but they have non-trivial quantum fluctuations. Accordingly, the beta func- 
tions that govern the change of couplings under the RG flow are quantum operator equations not 
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FIG. 4: (a) In the conventional RG, multi-trace operators are generated at low energies although only 
single-trace operators are turned on at UV. Once the initial condition is given, there is a unique RG trajectory 
determined by the classical beta function, (b) In the holographic description, one only needs to keep track of 
single-trace operators under the RG flow at the expense of making the sources for the single-trace operators 
dynamical variables. The partition function is given by sum over all possible RG trajectories for the single- 
trace operators. One has the freedom to employ different RG schemes, namely, one can choose different 
'speed' of RG flow at different scales. This freedom corresponds to the diffeomorphism invariance in the 
bulk. 



classical equations, 



(101) 



where we have chosen the gauge a (2;) = 1. The equation for the sources has to be supplemented by 
the equation for the conjugate operators, -j^^^ = [Pn, "H]. In the large N limit, the quantum beta 
function reduces to the classical equation in Eq. (|74] ) with the identification [A, B] -^i^, B} 
where {A, B} is the Poisson bracket. The 'Hamiltonian' that governs the quantum RG flow is 
dynamically generated from integrating out high energy modes at each step of RG. This is in 
contrast to the conventional RG where the RG trajectory is uniquely determined once an initial 
condition is given. In conventional RG, there is no quantum fluctuations for coupling constants, 
but one has to keep all multi-trace operators along the RG flow. This is illustrated in Fig. |4] 
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X. SUMMARY AND DISCUSSIONS 



From a first-principle construction, it is shown that a D-dimensional matrix field theory is 
mapped into a(D + 1) -dimensional quantum theory of gravity, where the metric in the bulk space- 
time is fully dynamical. The construction starts from the observation that one can identify high 
energy modes as fluctuating sources for the low energy modes in RG||8|]. For matrix field theories, 
this is implemented by introducing a dynamical source and its conjugate momentum for each pri- 
mary single-trace operator to remove multi-trace operators at each step of RG lllql . In particular, 
there is a spin two source and its conjugate momentum that represent the dynamical metric and 
the energy-momentum tensor, respectively. While the dynamical sources and momenta are ini- 
tially introduced as auxiliary fields, they acquire non-trivial dynamics as high energy modes are 
integrated out. On the one hand, the double-trace operators that are generated from single-trace 
operators through quantum correction provides the quadratic kinetic term for the conjugate mo- 
menta. On the other hand, the potential terms, including the curvature term for the /^-dimensional 
dynamical metric, are generated from the source dependent determinant for the high energy mode 
that is integrated out at each step of RG. The kinetic and potential terms together can be viewed 
as a (Z^ + 1) -dimensional action written in the canonical formalism, once the extra dimension 
corresponding to the length scale of RG is interpreted as a time. The bulk theory takes the form 
of quantum theory of gravity coupled with matter fields of various spins. Because of the freedom 
to choose different local RG scheme without modifying the partition function, one has {D + 1)- 
dimensional diffeomorphism in the bulk, which in turn leads to (D + 1) local constraints. The 
Hamiltonian constraint originates from the gauge freedom in choosing the spacetime dependent 
speed of coarse graining in the local RG procedure, while the D momentum constraints are asso- 
ciated with relabeling the D-dimensional coordinates of low energy modes relative to the coordi- 
nates of the high energy modes. Because different choices of local RG scheme merely correspond 
to choosing different gauge, the (D + 1) local constraints are first-class. 

The holographic dual for the matrix model includes dynamical gravity and other fields. Gener- 
ically, the cosmological constant in the bulk is expected to be order of the UV cut-off of the 
boundary field theory. Then the saddle point geometry in the bulk will have a curvature that is 
comparable to the scale that controls the gradient expansion for the action in the bulk. In this case, 
there will be no sense of locality within the distance scale over which the bulk spacetime is flat, 
although the geometry is classical due to the suppressed quantum fluctuations for a sufficiently 
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large A^. It would be of great interest to find boundary field theories whose gravity dual have a 
weakly curved bulk spacetime through the explicit construction. This would require stabilizing 
the theory at a strong coupling. 

For a general field theory, it is not easy to derive the dual theory in a closed form because one 
has to keep a large number of fields in the bulk. However, we have a concrete prescription to 
identify gravitational duals starting from boundary field theories. Using this prescription, one can 
try to examine the properties of the field theories which have simple gravity duals. For example, 
it will be interesting to see if one can identify the field theory whose holographic dual is the pure 
gravity. 
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XII. APPENDIX A : EXISTENCE AND UNIQUENESS OF CANONICAL METRIC 

We constructively prove the statement that there is one and only one metric in which the 
quadratic kinetic term has the canonical form as in Eq. (fTTT i for a given set of sources. Under 
a change of the metric used in covariant derivatives, only those operators that are quadratic in $ 
mix with the kinetic term. So we focus on the quadratic Lagrangian, 

oo 

^^'^ = -5Z/^'-^"tr[<l>9^,..9^„<l>]. (102) 

n=0 

The term with n = 1 can be absorbed into the term with n = via an integration by part, but it is 
more convenient to keep it for now. The goal is to find the metric in which the same Lagrangian is 
written in the canonical form, 

oo 

£(2) _ -J2^/\G\r"■■'^-tr[<l>V^,..V^„^, (103) 

n=0 
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where is the covariant derivative associated with the canonical metric that satisfies the condi- 
tion, G^^ = J^^ 

The canonical metric will be a local functional of the sources Our strategy is to 

compute the canonical metric using a gradient expansion of the sources. Suppose that G^^ is the 
metric that coincides with the canonical metric up to the v-ih. order in derivative. Namely, G^^'^ is 
made of the terms that have v or less derivatives of the sources in the canonical metric. The exact 
canonical metric is G^. The Lagrangian in Eq. (11021 ) can be expressed in terms of the operators 
constructed with the covariant derivative with the metric G^'^, 

oo 

To the zeroth order in derivative, the canonical metric is completely determined from j^" , 

V^Gr = 3'". (105) 

It is noted that Gq^ itself is completely fixed m D > 2. In this metric, the source for the two 
derivative operator becomes 

where Fq.^^ is the Christoffel symbol for the metric Gq'^, and ... represents the terms that have 
at least two derivatives of the source j^"^, such as j^°'^^^o-a'y^o-i3S j^°'^^^a^o-i3-y To the first 
order in derivative, the canonical metric is given by 

Note that the difference between G^'^ and Gq*^ has one derivative in the source. If we rewrite the 
Lagrangian using the covariant operators associated with Gi'^, the source for the two derivative 
operator tr[$V^V^$] differs from -/GiGi^ by terms that have at least two derivatives of the 
sources. If one repeats this procedure, one can uniquely determine G^*^ from G^^^^^ by adding 
terms that have v derivatives of the sources. This proves that there is a unique canonical metric for 
a given set of sources. 

There is a useful corollary. Suppose there are two theories which have canonical kinetic terms 
defined on two different curved backgrounds, 

Ca = -N^\G^ G'^^'tr [<I>V^V>] + (108) 
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with a = 1, 2. If Gi'^ ^ Gg'^, the two theories are distinct in the following sense. If the Lagrangians 
are re-expressed in terms of the operators defined on the flat manifold in Eq. (|2]), the two theories 
should have different sets of sources. 
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